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1. Introduction 


In this paper, we are interested in studying the following 2D 
nesq equations with fractional Laplacian dissipation 

" dtu + {u ■ V)u + + Vp = 6*62, X e M^, 

dte + {u-V)e + KA^e = Q, a:eM^ 

V ■ M = 0, X e 

^ u{x,0) = Uo{x), 6*(x, 0) = 6*o(x), x G M^, 


incompressible Boussi- 


f > 0, 
f > 0, 
t > 0, 


( 1 . 1 ) 


where the numbers i^>0,fi:>0,aG[0, 2] and /3 G [0, 2] are real parameters. Here 
u{x, t) = (ni(x, t), U 2 {x, t)) is a vector field denoting the velocity, 9 = 9(x, t) is a 
scalar function denoting the temperature, p is the scalar pressure and 62 = (0, 1). The 
fractional Laplacian operator A“, A := (—A) 2 denotes the Zygmund operator which is 
dehned through the Fourier transform, namely 


A“/(o = ier/(o, 


where 

(2vr)^ Jr2 

The fractional dissipation operator severs to model many physical phenomena (see CZI) 
in hydrodynamics and molecular biology such as anomalous diffusion in semiconductor 
growth (see [3S])- We remark the convention that by a = 0 we mean that there is no 
dissipation in fll.ip i . and similarly /3 = 0 represents that there is no dissipation in fll.lll o. 
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The standard Boussinesq equations (namely a = (3 = 2) are of relevance to study a 
number of models coming from atmospheric or oceanographic turbulence where rotation 
and stratihcation play an important role (see for example |32l |35]). Moreover, as point 
out in [32], the 2D inviscid Boussinesq equations, namely (ll.l|l with a = /? = 0 are 
identical to the incompressible axi-symmetric (away from the 2 :-axis) swirling 3D Euler 
equations. There are geophysical circumstances in which the Boussinesq equations with 
fractional Laplacian may arise. The effect of kinematic and thermal diffusion is atten¬ 
uated by the thinning of atmosphere. This anomalous attenuation can be modeled by 
using the space fractional Laplacian (see mm)- 

The global well-posedness of the 2D Boussinesq equations has recently drawn a lot 
of attention and many important results have been established. It is well-known that 
the system fll.ip with full Laplacian dissipation (namely, a = /3 = 2) is global well- 
posed, see, e.g., [7]. In the case of inviscid Boussinesq equations, the global regularity 
problem turns out to be extremely difficult and remains outstandingly open. Therefore, 
it is natural to consider the intermediate cases. Actually, many important progress has 
recently been made on this direction. Almost at the same time, Chae [10] and Hon and 
Li [23] proved the global regularity for the system fll.ip when a = 2 and /3 = 0 or a = 0 
and (3 = 2 independently. Since then, much efforts are devoted to the global regularity 
of fll.ip with the smallest possible a G (0,2) and (3 E (0,2). As pointed out in [21] . 
we can classify a and (3 into three categories: the subcritical case when a + (3 > the 
critical case when a + (3 = 1 and the supercritical case when a + [i < 1. As a rule of 
thumb, with current methods it seems impossible to obtain the global regularity for the 
2D Boussinesq equations with supercritical dissipation. Recently, Jin, Wu and Yang [25] 
established the eventual regularity of weak solutions of the system fll.ip when a and (3 
are in the suitable supercritical range. For the critical case, there are several works are 
available. In the two elegant papers, Hmidi, Keraani and Rousset [211122] established the 
global well-posedness result to the system fll.ip with two special critical cases, namely 
a = 1 and /? = 0 or a = 0 and (3 = 1. The more general critical case, that is a -1- = 1 

with 0 < a, /? < 1 is extremely difficult. The standard energy estimates do not yield the 
global bounds in any Sobolev spaces when a and (3 in the critical case. Very recently, 
the global regularity of the general critical case a + (3 = 1 with a > ~ 0.9132 

and 0 < < 1 was recently examined by Jin, Miao, Wu and Zhang [23]. This result 

was further improved by Stefanov and Wu [5B] by further enlarging the range of a with 
a + (3 = 1 and 1 > a > ~ 0.7981 and 0 < /3 < 1. Here we want to state 

that even in the subcritical ranges, namely a + (3 > 1 with 0 < a < 1 and 0 < /3 < 1, 
the global regularity of fll.ip is also dehnitely nontrivial and quite difficult. In fact, to 
the best of our knowledge there are only several works concerning the subcritical cases. 
More precisely, Miao and Xue [31] obtained the global regularity for system fll.ll) for the 
case z/ > 0, K > 0 and 


6-^6 


<«<!, 1 — a < (3 < min 


7 + 2\/6 ^ a{l-a) ^ 

a — 2, —;=-, 2 — 2a 


5 


VQ-2a' 


In addition, Constantin and Vicol [13] verihed the global regularity of the system fll.ip 
on the case when the thermal diffusion dominates, namely 


u > 0, K> 0, 0 < a <2, 0 < (3 <2, (3 > 


2^ a 
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Recently, Yang, Jin and Wu [15] proved the global well-posedness of the system fll.ip 
with 


> 0 , 


a 


K > 0, 0 < a < 1, 0 < /3 < 1, /?>!--, /3 > 


2 + a 
3 


/ 3 > 


10-5a 
10-4a' 


Very recently, the authors [H] established the global regularity for the 2D Boussinesq 
equations with a new range of fractional powers, namely z/ > 0, /t > 0 and 


0.783 


21 - yWf 


< a < 1, 


1 — a < /3 < min 


a (3a —2)(a; + 2) 2 — 2a'\ 

2’ 10-7a ’ 4a-3/' 


Here we also want to mention that the two works [T3l 145] have been improved by the re¬ 
cent manuscript mi- More precisely, the authors in |17] established the global regularity 
result for the 2D Boussinesq equations with 

Q/ 2 ~h 

z/>0, «;>0, 0<a<l, 0</5<l, /3>1 — —, /?> —-—. 

^ o 

The case of partial anisotropic dissipation has been considered in several settings (see 
for instance [I1E1IIS1I2S1I2S]). For the global smooth solutions to the damped Boussinesq 
equations with small initial datum, we refer the readers to the recent works [21 02]. 
Moreover, the global unique solution of the Boussinesq equations for the Yudovich type 
data has been established by many works, and we refer the readers to the interesting 
works It is worth remarking that there are several works concerning the 

global regularity for the 2D Boussinesq equations with logarithmical dissipation (see, 
e.g., |2ni 1121 EH] ) • Many other interesting recent results on the Boussinesq equations can 
be found, with no intention to be complete (see, e.g., [HI El El EH ESI ED 0310H] and 
the references therein). 

The goal of this paper is to establish the global regularity of solutions to the system 
fll.ip with a new range of fractional powers of the Laplacian. Since the concrete values 
of the constant z/, k play no role in our discussion, for this reason, we shall assume 
1 / = K = 1 throughout this paper. Now let us state our main result as follows 


Theorem 1.1. Let 0 < a < 1 and 0 < < 1 satisfy 


P* : = 


max < 
l3’ 

2 — a 
2 ’ 


4 — a^ 'I 
4 + 3a/’ 


2 

0 < a < -, 
2 

- < a < 1. 


( 1 . 2 ) 


Assume that {uo,9o) E x i7®(M^) for any s > 2 and satisfies 'V ■ uq = 0. Then 

the system admits a unique global solution such that for any given T > 0 


u E C'([0,T];77"(M2)) nL2([0,T];i7"+t(M2)), 


e E C'([0,T];i7"(M2)) nL2([o,T]; 77^+2 (r2)). 


Let us give some remarks about our result. 

Remark 1.2. On the one hand, one can easily check that > /?*, thus Theorem 
11.11 signihcantlv improves Theorem 6.1 of [13]. On the other hand, our theorem also 
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significantly improves Theorem 1.1 of [15], which obtained the global regularity result 
under the assumption 


(5 > max 


2 —a 2 +a 10 —5a 4 


Finally, we m have proved that the system fll.ip admits a unique global solution pro¬ 
vided 

a 2 — a 2-|-a 

13 > and /5 > — 

Obviously, Theorem 11.11 signihcantly improves the result of |17] . 


Remark 1.3. Through the proof, we hud that Theorem II.II is always true for (3 > 
with any 0 < a < 1. However, it is easy to check that when 0 < a < |, it holds 

(2 4 — a^4 2 — a 

l3’ 4 + 3a/ 2 

In fact, the proof of Case 1 (below) is much complicated than the proof of the case 
(3 > ^ with any 0 < a < 1. 


Remark 1.4. Through the proof of Theorem 11.11 we strongly believe that if one may 
establish Lemma ^27f \ under somewhat weaker conditions than (3 > (3*, then Theorem ll.il 
can also be improved. As suggested by Jin, Miao, Wu and Zhang in [2l], the expected 
subcritical result is/3>l — a(l — a< f3*). However, at the moment we are not able to 
weaken the conditions (3 > (3*. We will investigate this issue further in our future work. 


Remark 1.5. The nonlinear lower bounds for the fractional Laplacian [TB] or the Holder 
estimates for advection fractional-diffusion equations [3^ entails us that if one can show 
that for any given T > 0 

sup ||n(f)||cc, < c )0 or sup ||ci;(f)|| 2 < 00 , 

0<t<T 0<t<T 

under the assumption (3 > 1 — a, then the equations are well-posed in the smooth 
category up to time T. Here tc := V x u = dx-^U 2 — is the vorticity and C°‘ stands 
for the classical Holder space. 


We outline the main idea in the proof of this theorem. A large portion of the efforts are 
devoted to obtaining global a priori bounds for u and 6 on the interval [0, T]. According 
to the dehnition of (3*, the proof of Theorem 11.11 is divided into two cases, that is, 

2 r 2 4 - a^ 4 

Case 1 : 0 < a < -, max < -I < /3 < 1, 

- 3’ I 3 ’ 4 + 3a/ ^ 

2 2 — a 

Case 2 ; - < a < 1, -< d < 1. 

3 “ 2 

To start, let us say some words about the proof of the work 135], where the main idea 
of the work 051 is to consider the combined quantity G (see (I2.10h for more details) 

dtG + {u ■ V)G + A“G = -[7^^, u-V]9 + (1.3) 

Here and in sequel, we have used the standard commutator notation 

[JZ/s, u ■ V]0 := Tljiiu ■ '\/9) — u ■ '\/TZj39. 
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Invoking some commutator estimates and some computations, the combined quantity G 
satisfies 

sup ||G'(t)||i2+ [ ||A^G(r)||^2dr < cx), 

0<t<T Jo 

which is true for fj > ^ and fJ > 

Then they show the estimate 

sup ||G(t)||LPO < CO 
0<t<T 

for some po > 2. This estimate together with the iterative process leads to 

sup ||G(t)||LP < CO, for any po < P < oo, 

0<t<T 

which is valid for /3 > 1 - f, ^ ^ and ^ 

However, the main argument used here is completely different from the work [15]. For 
Case 1, in view of several commutator estimates, we can show by combining L^-norm 
of the combined quantity G and the temperature 9 

sup (||G||i, + ||A‘’0||iO(t)+ r (\\AiG\\h + \\A‘*h\\h)(T)dT<a=. (1,4) 

0<t<T Jo 

whenever 0 < p < |. 

For Case 2, by combining L^-norm of the vorticity uj and the temperature 9, one can 
conclude that 

sup (||w||i2 + ||A'^6'|||2)(f) + f (||A2a;||^2 + ||A^+2 0||^2)(r) dr < CO. (1.5) 
0 <t<T Jo 

whenever 0 < <5 < |. 

The above two bounds fll.4l) and fll.5l) are the key component of this paper. With the 
help of the two bounds (11.41) and (II.hh . we will establish the following key global bound 

max ||M(f)||Lr- < oo 

0<t<T 

for any 2 < r < co. Thanks to the nonlinear lower bounds for the fractional Laplacian 
established in [13], the following key estimate holds 

max IIV6l(f)||L°o < oo. 

0<t<T 

Finally, with the above estimate at our disposal, the global regularity of u and 9 following 
a standard approach (see for instance mmiEa). 

The rest of the paper is organized as follows. In Section 2, we obtain the a priori 
estimates for sufficiently smooth solutions of the system fll.ll) . Section 3 is devoted to 
the proof of Theorem 11.11 Finally, in the Appendix, we give the proof of Lemmas 12.11 
and 12.31 for the sake of completeness. 
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2. A PRIORI ESTIMATES 


This section is devoted to the a priori estimates which can be viewed as a preparation 
for the proof of Theorem 11.11 To simplify the notations, we shall nse the letter C to 
denote a generic constant which may vary from line to line. The dependence of C on 
other parameters is usnally clear from the context and we shall explicitly specify it 
whenever necessary. 

The first lemma concerns the following commutator estimate, which plays a key role 
in proving our main result. The proof can be performed by making use of the Littlewood- 
Paley technique. To facilitate the reader, we will sketch the proof in the Appendix. 

Lemma 2.1. Let ^ ^ ^ P ^ [2, oo) o-nd pi, p 2 G [2, oo]. Assume r G [1, oo], 

6 G (0,1), s G (0,1) such that s + (5 < 1, then it holds 

IIIAt/lsllflj, < C(p,r,i,s)(||V/||„||<,||B.«-. + ||/|U.||slliO' (2-1) 

In particular, 

IIIAI/Mbjs,. <C(p,r,«,s)(||V/||„||9llBS;'V + ll/lli^ll9lli”)- ( 2 - 2 ) 

Here and in what follows, Bp ,, stands for the classical Besov space (see appendix for its 
precise definition). 


In order to prove Case 2, we shall use the next two commutator estimates involving 
Up := a,,A-A 

Lemma 2.2 (see PI). Assume that ^ < ft < 1 and 1 < p 2 < oo, 1 < pi, p^ < oo with 
— + A _|_ A = X, Then for 0 < si < 1 — /5 and si -\- S 2 > 1 — ft, the following holds true 


F[}Zg,UG-V]edx < CIIA^i^lliPi ||A||w^2.p2||G'||lp3, 


(2.3) 


where uq '■= V'^A and denotes the standard Sobolev space. 


Lemma 2.3 (see [30] )• Tet ^ + ^ = p 2 < Pi, P 2 A oo and 2 < p < oo. Assume 

that 0 < ft < 2 and V ■ n = 0, then we have for any r G [1, oo] 

||[7^/3, r-V]0|Up <C(||Vr|Upi||0||51-;3 + ||r|Ub||0|U2). (2.4) 

P2.1 

For the sake of completeness, we will give the proof of Lemma 12.31 in the Appendix. 


Finally, let us recall the following fractional type Gagliardo-Nirenberg inequality 
which is due to Hajaiej-Molinet-Ozawa-Wang [TO] . 

Lemma 2.4. Let 0 < p, po, pi, q, qo, qi < oo, s, Sq, Si G M and 0 < -d < 1. Then the 
following fractional type Gagliardo-Nirenberg inequality 

holds for all v G Bp° ,j^ n if and only if 

71 71 71 

-s = (1 - d)(-So) + d(- si), s < (1 - d)so + dsi, 

p ''Po ''pi 

1 1 — d d 

- <-^-, if Po ^ Pi and s = (1 - d)so + dsi, 

q do qi 


(2.5) 
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1 1 - d d 

So 7 ^ Si or - < -1-, ij Po = Pi o,nd 

q Qo qi 

n n 1 1 — 'd -d 

So-T^s- or - < -1-, ij s 

Po p q qo qi 

A special consequence of H2.5\) is the following bound 


' ' 1 




il-A 


B: 


2,2 


'Sgo,oc(®2)’ 


where 2 — 2/5 < s < with ^ < f3 < 1. 


We also have 


s = (1 - 'd)so + {tsi, 

< (1 - '! 5 )so + iHsi. 

, 2/5-1 

2 - 2s’ 


( 2 . 6 ) 


^ ^ll'^"^ll3(R2)|b|liJ(l2), /5>0, 0<7<i (2.7) 

Remark 2.5. Lemma [2.41 is also true in the nonhomogeneous framework. 

It follows from the basic energy estimates that the corresponding solution (n, 9) of 
the system fll.ll) obeys the following global bounds. 

Lemma 2.6. Assume {uo, 9o) satisfies the assumptions stated in Theorem, li.il Then 
the corresponding solution {u, 6) of admits the following bounds for any t > 0 

mt)\\l^+ f\\A^0{T)\\l.dT<\\eo\\h, 

Jo 

\\0{t)\\LP < \\0 o\\lp, Vpe[2,CX)], 

+ f dr < (IlnolU^ + tPoh^r. (2.8) 

Jo 


2.1. Case 1. As well-known, when 0 < a, /5 < 1, it is impossible to obtain the global 
bound of (u, 9) by direct energy estimate method. Actually, applying operator curl 
to the hrst equation in fll.ip . we have the following vorticity w = V x u 

dtw + {u ■'V)w + A°'w = dxi9. (2.9) 

However, the ’’vortex stretching” term dxj^9 appears to prevent us from proving any 
global bound for w. To circumvent this difficulty, a natural idea would be to eliminate 
dxj^9 from the vorticity equation. To this end, we generalize the idea of Hmidi, Keraani 
and Rousset [2I1122] to introduce a new quantity. More precisely, we set the combined 
quantity 

G = u:-ng9, ng-.= dx,A-^, 

which obeys the following equation 

dtG + {u ■ V)G + A^G = -[JZg, u-V]9 + (2.10) 

Since u is determined by u through the Biot-Savart law, we have 

u = V^A-^oj = V^A-\G + ng9) = V^A-^G + V^A-^ng9 := ug + ue. (2.11) 

We are now in the position to derive the following estimates concerning G and 9. 
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Lemma 2.7. Assume {uq, 9q) satisfies the assumptions stated in Theorem \1.1[ Let 
{u,9) be the corresponding solution of the system U.l\) . If (3 > max||, then the 

following estimate holds for any max | < g < ^ and t E [0, T] 

\\G{t)\\l, + Wmrnl, + f (IIA^GIIi. + ||A^+^0||iO W dr < C{T, uo, 9o), (2.12) 

Jo 

where C(T, uq, 9q) is a constant depending on T and the initial data. 


Remark 2.8. Although the above estimate (12.121) holds for max < g < 

yet by energy estimate (12.81) and the classical interpolation, we hnd that fl2.12p is true 
for any 0 < p < |. 


Proof of Lemma \ 2.7\ Applying A^ to fll.ip o and taking the inner product with A^9, 
we obtain 




A^{u-\/9)A‘^9dx. 


(2.13) 


Hence, an application of the divergence-free condition and commutator estimate (12.21) 
directly yields 


A‘^{u-V9)A<^9 dx 


[A^,w V]0 A‘^9dx 


I V ■[Afu]9 A‘^9dx 

1R2 

< C||A'-«|A«,t,)9|U,||A«+««||„ 

< C|||A',«|9|| . ,||A«+I«||i. 


D- ■ 

-^2,2 


< \\u\\l2\\9\\l2)\\A^+^9\\l2 


B 


< C'(||a;|U2||0|Uoo + ||n|Udl^lUOI|A^+ = 0||L2 

— ^(ll^lU^ll^lli^oo -|- 1177^011 ll^ll^^oo ||n||L2| 


< ^\\A^^^9\\l. + C{l + \\G\\l2 + 




where we have used the following facts 


"<2 


Il=)I|a^+=0||l^ 

(2.14) 




and \\n09\\L2 <G\\9\\^p, ^- 3 ' 


Inserting the above estimate in (12.131) . we thus obtain 

14||A'>«||i. + < C(1 + IIGIIi, + iwr^p. 


(2,15) 






























9 


In order to close the above inequality, we need to consider the equation fl2.10p . To this 
end, we multiply the equation fl2.10p by G to obtain 




Gdx 


Gdx 


[JZp, u ■ V]6* G dx 
[JZ( 3 , ug ■ V]6* G dx 


— / [JZp, ue ■ V]0 Gdx. 


(2.16) 


Bounding the hrst term at the R-H-S of (I2.16p according to the Holder inequality and 
the interpolation inequality, we thus get 

f A^-^d,,9Gdx < G\\A^+'^-^9 \\l2\\A^G\\l2 

< C||e||iri|A'+f9||I.||A*G||„ 


2 + q; — 3/3 2 + CK — 2/3 . , 

g > -- T = -- e (0, 1) 

^ 2 2p + /3 ^ ^ 

< -||A2G||^2 + q\\A^~^^9\\12 + C'||6'||i2. 


(2.17) 


Next we appeal to the commutator estimate fl2.4p and the fractional type Gagliardo- 
Nirenberg inequality fl2.6p to bound the third term at the R-H-S of fl2.16p 

[ [JZf,,ue-V]9Gdx < \\[n^,ue-V]9\\L2\\G\\L2 

Jr'2 

< GiWVuehAlOhi-^ + \\uo\\^^\\9U2)\\Gh2 

< q|A^-/^0|U4||0||^i-.||G|U2 + ||0||i2||G|U2 


4,1 


< C’|| 9 ||'A-,||Gn, + || 9 „||i,||Gn,, 


(2.18) 


where we have used the classical embedding Thanks to fl2.6l) . we have 

that for 2 — 2/3 < s < with ^ < /3 < 1 




A = 


2/3-1 
2 - 2s 


e (0, 1). 


According to Sobolev interpolation, we can get that for any ^ < s < ^ + g 


ll^ll 


B. 


2,2 


< 


|l-z 


^-7 


I|A^^ = ^IIl2(m 2), / = ^e(o, 1). 


Q 


(2.19) 


( 2 . 20 ) 


Inserting fl2.20p into fl2.19p and considering fl2.18p . we can conclude 

JjV.f,u,-V]eGdx < Glidin'.^,,,,||0|fy7,yi|G|u. 

< riimi|2A(l-Z)|| A^l+4nll2A^| 


C||«||“^‘-'’||A«+Sfl||B'||»„fy7,">,||G||„ 

1 B 2X(l-l) 

< «l|A'’+i«|li.+C||(l|| V" (l + ||G||i,), 

o 


( 2 . 21 ) 
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where in the last line we have used the following fact 


0 < s < 


f (2/3 — 1) + p 


2/3 - 1 s - f 


A/ < — 

213-1 + Q “2 2-2s g 


Combing all the restrictions on s yields 

/3' 


max <0,2 — 2/3, 

I z 

which would work as long as 


< s < mm 


< - ^ - < 2 . 

-2 l-Xl - 


f (2/3 — 1) + £) 


3-2/3 ^ ^_ 

2 ’ 2 2p-l + g 


Q > 


4-5/3 




Now we focus on the second term at the R-H-S of fl2.16p . The estimate fl2.7p as well as 
energy estimate ( 12 .Sp leads to 


l|A"'’e|| 1 1 < c\+o\\^xm+l^ <<>=, 0 < 7 < 




( 2 . 22 ) 


The commutator estimate (12.31) with si = 7/3 < 1 — /3, 1 — /3 — 7/3 < S 2 < f i Pi = 


7 , P2 = 


— and p 3 = 


2+2s2—Q: 2+Q:—2s2—47 

I [TZg, ug ■ V]0 G dx 

Jm? 

< C||A"'’«|| .||G|| l|G|| LP3 


allows us to show 


' L^t 

< C'||A'’''«||^i||AfG|U,(||G||lVI|A»Gr 


L 2 ^ 


a - 47 / , 2s2 + 47 - a _ 

—-—- < S 2 < a — {a + 2)7 p = -;-G (0, 1) 


a 


< sl|AtG||i.+C||A’'’«||‘r||G||i. 


LT 

< i||AtG||i,+C(l + ||A'"’«||h)l|G||i„ 

o Lt 

where in the last line the following fact has been applied 

2 / 1 , 2 s 2 + 47 - a 

- < — ^ fi = -< 1 — 27. 

1-/17 O! 

Putting all the restrictions on S 2 together, we have 


(2.23) 


1 - /3 - 7/3 < S 2 < 


a 


a — 47 


2 ’ 2 

Consequently, the above S 2 would work as long as 


< S 2 < a — {a + 2 ) 7 . 


max < 0 


2-213-a 


which leads to the key assumption 


1 1-/3 a + /3-l- 

< 7 <mm<;-, —-—, --f, 

^ ' 2 ’ /3 ’ 2 + a-/3i’ 


/ 3 > 


4 — 0 ;^ 
4 + 3q; 
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It is worth noting that the fact > 1 — a and this is the only place where we use the 
assumption (3 > Inserting the above aforementioned estimates fl2.17p . fl2.2ip and 

(I2.23P into (I2.16p yields 

2^11^11^2 + ^11^^ G'lliz < -\\J\.‘^^^9\\\2 + C'll^ll^z + (7(1 + ||0||^|)(1 + ||(7||^2) 

+C{l + \\A^f^ef\)\\G\\l2. (2.24) 

Summing up fl2.24l) and fl2.15p . we thereby obtain 

4(||G||i. + ||A«»||i.) + ||A»G||i. + ||A«+«9||i. < (7(1 + ||9||% + ||A’'’«||7)(1 + l|G||i,). 

at L~i 

which together with the classical Gronwall inequality and (I2.22p lead to 

sup (||G(()||i. + ||A«»(()||i,) + r (||A«G||i, + \\A‘*h\\h)(T)dT < C(T, u„, «„). 

0<t<T Jo 

This completes the proof of Lemma I2.7[ □ 


2.2. Case 2. In this case, we consider the vorticity uj instead of the combined quantity 
G. Now we derive the following estimates concerning vorticity u and the temperature 
6, which can be stated as follows. 


Lemma 2.9. Assume (uq, 9o) satisfies the assumptions stated in Theorem M . R Let {u, 9) 
he the corresponding solution of the system U.l]) . If fi > then the vorticity u and 
the temperature 9 admit the following bound for any 5 < | and t e [0,T] 

II^(^)IIl 2 + + [ (II^^^IIl 2 + ||A‘^+20||^2 )(t) dr < (7(r, n-o, ^o)i (2.25) 

Jo 

where C(T, uq, 9o) is a constant depending on T and the initial data. 


Remark 2.10. Similarly, by energy estimate fl2.8p and the classical interpolation, we 
hnd that (12.251) is true for any 0 < h < |. 

Proof of Lemma \2.iA With the same argument used in obtaining fl2.14p . we hnd that 


\fL 

2 dt 


\\kH\\l. + \\k‘*^e\\l. < G||A‘-f[A^u]0||„||A'+i9n, 

< G|||A',u|9||^,_,||A'+?9n. 

< G|||A‘,t.|9|| ._j||A‘+l9|U. 

^2,2 


< c(iivun.|i«ii + iiu||g|i«iii>)iia*+>«iu. (i < I) 

00,2 

< G{\\u\\lA\0\\l^ + \\u\\lA\0\W)\\I^'^-^9\\l2 

^ d||^2 + (7||do||ioo||w||^2 + (7||'u||i2||d||^2. 
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Substituting the above estimate into fl2.26p . we thus obtain 




Mh- 


(2.26) 

In order to obtain the global bound of the velocity u, we resort to the vorticity w 
equation fl2.9p 

dtw + {u ■ V)w + A°'w = dx^O. (2.27) 

Testing it by u yields 


1 d 

2 dt 


Icnlli. + ||A2n;||^, = 


dxi 9 u dx 


< C||A'-t0|U2||Afa;|U2 

< C7||A^0||i7-||A^+^0||£.||Afn;|U2 

< -|| A 2 a ;||^2 + 6\\\2 + C || A ^ 0 ||^ 2 , 


(2.28) 


where we have applied the following Sobolev interpolation 


||A^-f 0 IU 2 < C||A^0||ir||A^+^0||E2, 


2 - « - 2h 


r = 


Note the fact 


a-/3 , 2 

■ < 6 < 


a 


2 ' ' ' 2 

Putting all the restrictions on b together, we have 


0 < r < 1. 


€ ( 0 , 1 ). 


2 — a — fb ^ (2 — a (3 

- < 5 < mm <-, — 

2 I 2 ’ 2 


2 ' 


Thus, this is the only place in the proof where we use the assumption of the theorem, 
namely j3 > 

Summing up (12.281) and (I2.26p . we thereby obtain 

^(ll^llia + ||A^0||^2) + ||A2a;||^2 + HA^+a^H^^ < C(||a;||^2 + ||A^6'||^2) + C'||'»-Ili2||0||^2, 

which together with the classical Gronwall inequality leads to 

sup (||n;||i 2 + ||A^ 0 || 22 )(t) + / (||A^a;||i 2 + ||A‘5+20||2^) (r) dr < cx). 

0<t<T Jo 

This completes the proof of Lemma 12.91 □ 


Both in Case 1 and Case 2, we can establish the following global a priori bound 
u(t)\\Lr for any 2 < r < 00 and 0 < t < T. 


Lemma 2.11. Assume 13 satisfies the assumptions stated in Lemmas \2. 7| and \2.tA then 
the velocity field u obeys the following key global a priori bound for any 2 < r < 00 and 
0<t<T 


sup \\u(t)\\Lr < C{r, T, Mo, ^0), 

0<t<T 


where C{r, T, uq, 9o) is a constant depending on r, T and the initial data. 


(2.29) 
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Proof of Lemma \2.11[ Let us notice that in Case 1, we have /^ > §• As a result, we 
can select q satisfying 1 — (3 < q < ^ such that 

WWl^ < ||A'-^ 0 ||l 2 < \\6\\h. < oo. 

Recalling G = u — 7lp9 and the estimate 02.121) . we get 

||<^||l 2 < ||G||l 2 + ||7^^0||2,2 < oo, 

which together with 02.25P implies that we have both in Case 1 and Case 2 

sup ||a;(t)||i:,2 < oo. 

0<t<T 

By the Sobolev interpolation inequality 


sup IIn 

0<t<T 


i-i 


Lr < C{r)\\u\\l2\\Vu\\^^ 


< C(r)||n||£2||a;||].2" 


< C{r, T, uo, 6o). 

Consequently, this immediately completes the proof of Lemma 12.111 


(2.30) 

□ 


With the help of the above estimate 02.29p . we are able to the next lemma, which is 
concerned with the global a priori bounda || V0||i,oo as well as ||a;||i,oo. 


Lemma 2.12. Assume {uo, Oq) satisfies the assumptions stated in Theorem M . 1\ Assume 
{3 satisfies the assumptions stated in Lemmas \2. 7| a,nd \2.!A then the temperature 9 and 
the vorticity u admit the following key global a priori bound 


sup ||V6'(f)||L<x. < C{T, Uo, 9o), 

0<t<T 

(2.31) 

sup ||a;(f)||L<x. < C(T, uq, 9o), 

0<t<T 

(2.32) 


where C{T, uq, 9q) is a constant depending on T and the initial data. 


Proof of Lemma 12.121 The idea of the proof is based on the argument of nonlinear 
lower bounds for the fractional Laplacian established in [I3] . For convenience the reader, 
we present the details as follows. By the elementary calculations, it is not hard to check 
that 

fi* > ^ 


1 -|- 


Therefore, this fact further implies 




1 T cr 


We start with the following pointwise bound 


V/(x) • A“V/(x) > -A“(|V/(x)r) 


ap 

P+2 

LP 


(2.33) 
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which can be proved by combining the proofs of Theorems 2.2 and 2.5 of |13j . 

Applying V to the temperature equation of fll.ip and multiplying the resulting equation 
by lead to 

]^{dt + u ■ v)|V0p + ve ■ kNe = -Vu ve ■ ve. (2.34) 

Thus, by making use of fl2.33p with p = oo, we immediately arrive at 

l-{dt + u-V + A^)\ve\^ + ^ ^ 235 ) 

Suppose that |V6*(x,f)| achieves the maximum at the point x = x{t), then we get 

a,|V0(i,i)P + c,2k^ < 4>(()2 ||Vu||i«. (2,36) 

ItoIIloo 

where 

$(f) = ||V0(.,f)|U.. 

Similarly, let us assume that |a;(x,f)| achieves the maximum at the point x = x{t) and 
denote 

Q{t) = ||a;(.,t)||Loc. 

Recalling the vorticity equation 

dtw + {u ■ V)w + A"zc = 

and adapting the same argument used above, we can conclude that 

dt\uj{x,t)\^ + C 2 -< $(f)f2(t), (2.37) 

where the number r G (2, cxo) will be fixed hereafter. 

To bound we need the following logarithmic inequality which was established 

in ([13]) 

II Vm(., f)||Loo < Co + C'of2(t) 

+Con(f) log++ J (l + A'(r) + n(r) + <h(r))’" drj (2.38) 

where Co > 0 is a constant depending on initial data, K{t) is a bounded function on 
the interval [0, T] and T = T{a,(3). 

Therefore, it follows from fl2.36p and fl2.37p that 

dt\we{x,t)\^+ cMtf+^ < Co<h(f)2|i + f2(f) + fi(f)iog+(i 

+ ^ (l + i7(r)+f2(r)+ <h(r))^dr)}, (2.39) 

dt\uj{x,t)\^ + (2.40) 

where <h(f) = \V6(x,t)\ and Q{t) = |a;(x,f)|, and the constants Co, Ci, C 2 depend on 
the initial data, a, j3 and ||M(T)||j,r. 

Suppose that M > 0 is large enough to be hxed hereafter. Assuming the solutions blow 
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up at time T, thus = cx). Then we can select Tq G (0, T) as the hrst time 

such that <h(To) = M > 4<h(0). Now one may deduce from fl2.40p that for any t G [0, Tq] 


Q{t) < max 


(f) 


2 + r 

2 + (l + Qi)r 


/ M \ 2+{l'+a)r- 




:= M, 


as long as M is large enough in terms of f2(0), a, r and C 2 . Let us give details about 
how to get the above estimate. Actually, if Q{t) > M, then 

> C'2fi(t)2+^ - MQ{t) 

> - M)Q{t) 

> (C2M^+^ - M)Q{t) 

= 0 . 

Thus it follows from fl2.40l) that dt\oj{x,t)\'^ < 0. This implies that fl{t) cannot exceed 
the value M. Hence, the following inequality is an easy consequence of (I2.39h 

dt\ve{x,t)\^ + < C'o$(t)2|l + M + Mlog+ {l + {l + K{T) + M + Mf^ 

Repeating the same argument as above, we obtain 


< 


max 


^(l + M + Mlog+(^l + (l + iL(r) + M + M)’ 


2 + r 


for any t G [ 0 , Tq] 

Now notice that M ^ M 2 +(i+“)’-, we select M large enough such that 

C*! 

which is equivalent to 


l + M + Mlog+ ( 1 + {I + K{T) + M + 


2 + r Mf^ 

X _|_ ]\^2 + {l + o,)r (^X + log_,_ M) < 

(_y 


(2.41) 


Thanks to the fact {3 > -p^, it is sufficient to choose r as 

2 ( 1 -/?) 


1+0’ 
ro < r < 00, 


Tq = max 


2 


(1 + cr)/? — 1 

then the above inequality (I2.4ip can be guaranteed due to the following fact 

2 + r 


ro < r ^ 


</?. 


2 + (1 + a)r 

Hence, it is not difficult to verify that $(To) ^ X’ ''^^ich contradict the dehnition of 
Tq. We thus get the fact that $(t) never blows up as t ^ T when T < 00. As a direct 
consequence of above fact, we infer that 

sup ||V 6 '(t)||Loc < C{T, uq , 6 >o) < 00. 

0<t<T 


As a consequence of the above estimate, it follows from the vorticity equation fl2.27p 
that for any 0 < t < T 

\\uj{t)\\L<=- < llc^olU- + [ ||V6'(r)||ioo dr < C{T, uq, Oq) < 00 . 

Jo 
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This concludes the proof of Lemma I2.121 


□ 


3. The proof of Theorem 11.11 


In this section we give the proof of Theorem ll.il With the estimates fl2.3ip and fl2.32p 
at hand, the proof can be performed as the classical approach. 


Proof of Theorem \1.1{ To begin with, we smooth the initial data to consider the 
following approximate system 

x e t > 0, 

=0, xe 


V • = 0, 

^ = SnUo{x), = SnOo{x), 


X G 
X e ] 


P, f > 0, 
t > 0, 


(3.1) 


where S'at is the low-frequency cut-off operator (see Appendix for it dehnition). 

Now we apply (/ -|- A)® to system fl3.ip and multiply the resulting equations by 
(/ -|- AYu^^'^ and (/ -|- AyO^^'^ respectively, add them up to to conclude that 

< C (1 + iivm^iU" + 

< C(l+||«l''l|l«- + ll«<"'ll«-)(iy''’llff. + ll«<"'ll?,.). (3.2) 

where we have used the embedding H^(E?)) for any s > 2. 

Therefore, there exists a time 

T*:=C-(l + (||«„||?,. + ||e„||?,.))"" 

for some absolute constant C* > 0 such that 


^{N) ^ l°°{[o,t*y hyrY) n lY[o,t*y h^+^{rY), 

q{N) ^ ^ooQg^ p 

Note that 

- A“n(^) + V9^^'>e2, 

where V denote the Leray projection onto divergence-free vector fields. 
Thus, it is not hard to see that 

G L~([0,T*); H^-\rY). 


Consequently, we assume that 

at.!"', a,#!"’ e iLdo.r-); ^^r'(R")). 

Since the embedding is locally compact, the well-known Aubin-Lions argu¬ 

ment and Cantor’s diagonal process, we conclude that there exists a solution satisfying 

















17 


e e L°°([0,T*);//*(M2)) nL2([0,T*);i/"+^(M2)). 

The continuity of u and 9 in time, namely u, 6 E C{[0, T*); can be obtained by 

a standard approach. It suffices to consider u G (^([O, T*); as the same fashion 

can be applied to 6 to obtain the desired result. 

By the equivalent norm, it yields 

ll«(ti) - = {(53 + 53)(y‘||A,«((,) - A,«(t2)|U,)2}f (3.3) 

j<N j>N 

where Aj is the non-homogeneous Littlewood-Paley operator (see Appendix for its def¬ 
inition). Let £ > 0 be arbitrarily small. Due to n G L°°([0, T*); there exists a 

integer > 0 such that 

{ 53(2'‘||A,u(ii) - A^«(i2)n.)2}’ < (3.4) 

j>N 


Recalling the system fll.ip i . we obtain for 0 < ti < t 2 < T* that 


Aju{ti) - Aju{t 2 ) 



A. 

dr 


Ajuir) dr 



AjV[ 6 e 2 — {u ■ V)m 


A"n](r) dr. 


Therefore, we can get 


(3.5) 




Aj'it(fi) - Aju{t 

j<N 

= - (w V)n-A"n](r)dr 

j<N 

< r|||A/|U. + ||A,(«.V«)|U. + ||||A,A“«||y](r)<ir 

j<N “^*1 

= r + 2^(^-')||A,(n ■ Vu\\l 2 + 2^(*-i+“)|| A,n|U 2 ](r) dr 

j<N 

< C y~^ 2^-^ ^||0||^s-i |ti — ^ 2 ! + II (11 ■ V)n||^s-i |ti — ^ 2 ! + ||m||^s-i+c. |fi — ^ 2 !^ 

j<N 

— ^ 'y ^ — ^ 2 ! ^II^IIh«-i a II^IIl°° II^^II^'’-! a II^^IIl°° II^IIh®-! a IIi^IIh®^ 

j<N 

< (A22^|fi-f2|(||ii|iy A IIwIIh^ + ||0||h^), (3.6) 


where the Sobolev imbeddings iL®“^(M^) and iL®“^(M^) L°°(M^) with 

s > 2 are used several times in the last inequality. 

Thus, the following holds true 

{ 53(y‘||A^u((i) - A,«((2)|U,)2}’ < I (3.7) 

j<N 
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provided \ti — t 2 \ small enough. 

Combining fl3.4p with fl3.7p implies u G C{[0,T*)] Moreover, the uniqueness is 

clear since the velocity and the temperature are both in Lipschitz spaces. 

Now, it remains for us to show that the local smooth solutions may be extended to 
all positive time. It suffices to state that under the assumption of the theorem and any 
given T > 0, we have 

sup (llu(i)lll,. + ll«(()ll?,.) + r + ll«(*)ll" ») * < C{T, »„). 

0<t<T Jo H 

In consequence, the energy estimate (13. 2 p ensures that 

< c(i + iivuiu» + iiv»iu.)(i|u|iy + ii»iiy). 

To obtain the global existence of smooth solutions, the standard procedure is to bound 
the term ||Vn||Loo with ||cii||l°o and a Sobolev extrapolation inequality with logarithmic 
correction (see e.g., [HIE]) 

II Vn||Loo(]R 2 ) < c{l + ||tt|| 2 , 2 (]R 2 ) + ||a;||ioo(]R 2 ) ln(e + ||ti||/^s(]K 2 ))^, s > 2. 
Consequently, it enables us to get 

|(ii«wiiy+ii9(<)iiy) + 

— C{\ + ||a;||Loo + II V0 ||l°°) ffi(e + ||ii||ir® + ll^ll^s)(||'^llrfs + ll^llri®)- (3-8) 

Applying the log-Gronwall type inequality as well as the estimates fl2.3ip and (12.321) . we 
eventually obtain the desired estimates. This concludes the proof of Theorem 11.11 □ 


Appendix A. The proof of Lemmas 12.11 and 12.31 

Before proving Lemmas 12.11 and 12.31 we hrst recall the so-called Littlewood-Paley 
operators and their elementary properties which allow us to dehne the Besov spaces (see 
for example [3l H] [33l [39|). It will be also convenient to introduce some function spaces 
and review some well-known facts. 

Let (x, (/?) be a couple of smooth functions with values in [0,1] such that x G C((°(M"') 
is supported in the ball S ;= {^ G M", |^| < |}, G Cg°(M'^) is supported in the annulus 
C G M", f < 1^1 < §} and satisfy 

+ = VfeB", ^45(2-y) = l, V5eR"\{0}. 

isN jez 

For every u G S" (tempered distributions) we dehne the non-homogeneous Littlewood- 
Paley operators as follows, 

AjU = 0 j < —2; A_in = x(-D)m; Vj G N, Ajtt = ^p{2~^D)u. 

We shall also denote 

SjU := AfcU, AjU := Aj_in -|- A^n -|- Aj+iu. 
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We now point out several simple facts concerning the operators Aj: By compactness of 
the supports of the series of Fourier transform, we have 

AjAiu = 0, \j — l\>2 and Ak{SiuAiv) = 0 |A: — /| > 5. 

for any u and v. Moreover, it is easy to check that 

snppT{Sj-iuAjv) ^1^1 ^2^ < 1^1 < 
supp I 1^1 < 8 X 2-^1 , 

where T denotes the Fourier transform and B to denote C~^B < A < CB for some 
positive constant C. 

Let us recall Let us recall the dehnition of homogeneous and inhomogeneous Besov spaces 
through the dyadic decomposition. 

Definition A.l. Let s G M, (p, r) G [1,+cxo]^. The inhomogeneous Besov space B*,, are 
defined as a space of f E S' (ML) such that 

By = {/G5'(M"),||/||B^,, <oo}, 

where 

f (5^2'"||A,/r„)t Vr<oo, 

II/IIb- = < e-1 

I sup 2^^||Ay||iP, V r = CX). 
f i>-i 


Many frequently used function spaces are special cases of Besov spaces. For s G 
M, (p, r) G [1, +C)o]^, we have the following fact 


For any s G M and 1 < g < oo. 



\H’> 


^g,min{g,2} 




^ B 


S 

g,max{g,2} ‘ 


Bernstein inequalities are fundamental in the analysis involving Besov spaces and 
these inequalities trade integrability for derivatives. 


Lemma A. 2 (see 0). Let k E N U {0}, 1 < a < fo < oo. Assume k = \a\, then there 
exist positive constants Ci and C 2 independent of j and f only such that 

suppfc {{ e E” : |{| < V} ^ || 8 "/|U. < Cl 2 i‘+C(i-i)||/||i,.; 

snppfc {? € E” : ICI SI 2'} ^ Ci2'‘||/||i. < ||a"/IU» < C2 2J‘+«J-t)||/||i.. 

Here we use A < B to denote A < CB for some positive constant C. 


To prove Lemmas 12.11 and 12.31 the following lemma will be used extensively. 

Proposition A.3. Given (pi, P2) G [2, 00]^ and p E [2, 00) such that ^ ^ 

f, g and h he three functions such that V/ G L^L 9 ^ o,f^d xh E L^. Then it holds 

\\h-k{fg) - f{h-kg)\\LP < ||x/i||ii||V/||LPi||^||LP 2 , 
where x stands for the convolution symbol. 
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Proof of Proposition IA.3i We remark that Proposition lA. 31 with pi = p and p 2 = oo 
has been proven in [22]. The interested reader may refer to [23] for general case. To 
facilitate the reader, we give the details. 

By direct calculation, one may easily show that 


h-k{fg){x)-f{h-kg){x) = / h{x-y)g{y){f{y)-f{x))dy 



- y)9{y){y - a^)-(V/)(x + {y- x)t) dydt 

- |)^-(V/)(x - z)^dzdt. 
According to the Minkowski inequality and the Holder inequality, one has 



\\h-k{fg) - f{h-k g)\\LP < 

< 

which is nothing but the desired result. 



^(D^IIWIIlpII^IIlS^ dzdt 


E2 Jo 

< ||x/l||ii|| V/||lpi ||fi'||LP2, 


(A.l) 

□ 


Now let us proceed to prove Lemma [2T1 To start, we use Bony’s decomposition to 
present the commutator as 

^k[A\f]g = A,([A^ Vi/]^.^)+ E Afc([A^ A,./]Vw) 

|j-fc|<4 It-fc|<4 

+ A»(|A', A,/]A,j) 

j—k>—A 

:= Ni + N 2 + N^. (-^*2) 

Now we recall the following fact. Let A bet an annulus centered at the origin. Then 
for every F with spectrum supported on 2^A, there exists g G whose Fourier 

transform supported away from the origin, such that 

A^F = 2 A^+^^g{ 2 A)i,F. 

For hxed fc, the summation over \j — k\ <4 involves only a hnite number of j's. For 
the sake of brevity, we shall replace the summations by their representative term with 
j = k in Ni and N 2 . In view of the above fact, Berstein’s lemma and Proposition IA.31 
we thus get 

IliVilUp < C'|lx2'’("+^)r/(2'’x)|Ui||VFfc_i/|Upi II Afc5'||LP2 

< C2'=(^-')||V/|Upi||Afc^|Up2. (A.3) 


Similarly, one can also deduce that 


||A^2||lp < C||x2"(-+'5)r7(2"x)|Ui||AfcV/|| LPi ||'S'fc-lfi'||LP2 

< C2^(^-^)||V/|| LPi E ii^'pii LP2 

l<k-2 

< C\\Vf\\ LPI E 2('’-')(^-i)2'(^-i)||Az(7|Up2. 

l<k-2 


(A.4) 
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Finally, the last term can be rewritten as 

Ns = A,(a'(A,/A,-^)-A,/A'A,^ 

j—k>—4 

= Afc (a'(A,/A,^) - A,fA%g 

j-k>-A, j>0 

+ Y1 A,(a'(A_i/A_i^)-A_i/A'A_i^7 

-l-k>-4 
■ = Nsi + Ns2 


(A.5) 


By Berstein’s lemma, the term N^i can be bounded without using commutator structure 


IIAsiIUp < C' (||A,(A^(A,/A,r?)) 

j—k>—4, j>0 

< C Y. {\\4‘(^ifAs) 

j—k>—4, j>0 


LP 


+ 


Ak{AjfA^Ajg) 
A,fA%g 


LP 


LP 


< C Y. 2 >‘‘-«||AjV/|| LPl II Aj5'||iP2 

j—k'>—4^ ji'>0 

< c Y 2'<'’-‘'||V/|U„||A,9|h 


\LP2 ■ 


(A, 6 ) 


j—k>—4 

Resorting Berstein’s lemma again, the term A 32 admits the following bound 

IIA 32 IILP < II Afc(A'^(A_i/A_i5f)) IliP + II Afc(A_i/A'^A_i5') ll^^p 


(A.7) 


-l<fc<3 -1<A;<3 

< C Y (2“ + l)l|A_i/n,,||A_,9||i,. 

-l<fc<3 

A C'A{-l<fc<3}||/||L2|lfi'IU2. 

By the dehnition of we have 

||[A',/]^||b^, < ||2^*||Ai|Up||^, + 1|2^*||A2|Up||^. + ||2^^||A3|Up||,. 

p,r M II llt^ II lU^ 

< ^^l|V/|| LPl ||2fc(s+5-i)||A^^||^p^|| 

II lU^ 

+C'||V/||lpi I y 2 (''-')(^+^-') 2 '("+^-i)||Ai^|Up 2 (s + h - 1 < 0) 

I / J j^r 

l<k-2 

+C'||V/|Upi y 2('’-^>2A*+^-b||A,.^|Up, ^ (5>0) 

j-k>-A 

+ C^II/IU2|b||L2 

< (^(p,r,5, s)(||V/||lpi|| c/||^.+«-i + ||/||l2 11 ^ 11 ^ 2 ). (A.8) 

Therefore, this concludes the proof of Lemma [2.11 

The proof of Lemma 12.31 is very similar to that of Lemma 12.11 Indeed, we can regard 
the operator 77/3 as the operator A^~^ without any difference. Now we just view f = u, 
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g = V9 and S = 1 — /3. Therefore, we have 

Ak[n0,u-v]e = ^k{[ng, Sj.,u-v]A,9'^ + ^k{[ng, a,u-v]s,.i9'^ 

|j-fc|<4 Ii-fc|<4 

+ Ak(jTZg, AjU ■ \/]Aj9'^ 

j—k>—4 



:— Ni -|- N 2 -\- Ns- 

The same as Ni 

and N 2 , we can conclude 


IIYiIUp < C2'=(i-«||Vm|Upi||A,0|Up. 

and 



\\N2 \\lp<C\\Vu\\l., 2('-'=)^2'(i-«||A;0|| 


l<k-2 

However, we need to deal with the term differently. We rewrite as 

Ns = Y 1 A.-w V]A,0) 

j—k>—4 

= Ak(TZi3{AjU ■ \/Aj9) — AjU ■ \/'JZpAj9'^ 

j—k>—4, j>0 

— Afc ^7 ?.^(A_im ■ VA_i0) — A_iri ■ 

-l-A;>-4 

:= A 31 + Ns 2 - 


By Berstein’s lemma and the divergence-free condition 


IIA31IUP < c { 

j—k>—4, j>0 

Ak{TZp{AjU ■ VAj 6 ')) 

+ 

LP 

A^i^AjU ■ \/7^gAj9^ 

LP^ 

c 

W 

VI 

AkijZpV ■ (AjuAj9)) 

+ 

LP 

AfcV ■ [Aju7ljsAj9^ 

LP^ 


j—k>—4, j>0 


j—k>—4, j>0 

j—k>—4, j>0 


and 


||A'32||lp < C'X{_i<fc<3}||'u||Lr||6'||2,2. (A.9) 

Putting all the above estimates together and making use of the Young inequality for 
series convolution yield 

mg,u-V]9\\L. < ||A47^;5, V]0||lp 

k>-l 

< C\\Vu\\lpi 2"^'"''lAfc0|Up. 

k>-l 
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+C'||Vm|Upi ^ 

k>-l l<k-2 


+C*|| Vm||lpi E E 2(^-^')(^-^)2(i-/^)^’||Aj-6»||lp2 

k>-l j-k>-A, j>0 


+C||Vm|U„ 5^ 2'-'2<'-'=W||A/||« 

k>-l j-k>-A, j>0 

+C X{-l<fc<3}||'*^||L"-11^11x2 

k>-l 

< c\\Vu\\LPi\\e\\j,i-p + \\u\\LA\o\\L^. (A.io) 

P2.1 


This completes the proof of Lemma 12.31 
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